Energy partition at the collapse of spherical cavitation bubbles 
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Spherically collapsing cavitation bubbles produce a shock wave followed by a rebound bubble. 
Here we present a systematic investigation of the energy partition between the rebound and the 
shock. Highly spherical cavitation bubbles are produced in microgravity, which suppress the buoy- 
ant pressure gradient that otherwise deteriorates the sphericity of the bubbles. We measure the 
radius of the rebound bubble and estimate the shock energy as a function of the initial bubble 
radius (2-5.6 mm) and the liquid pressure (10-80 kPa). Those measurements uncover a systematic 
pressure dependence of the energy partition between rebound and shock. We demonstrate that 
these observations agree with a physical model relying on a first-order approximation of the liquid 
compressibility and an adiabatic treatment of the non-condensable gas inside the bubble. Using 
this model we find that the energy partition between rebound and shock is dictated by a single 
non-dimensional parameter ^ = Ap7^/[p^o^^^(pc^)^~"^^^], where Ap = poo — Pv is the driving pres- 
sure, Poo is the static pressure in the liquid, pv is the vapor pressure, pgo is the pressure of the 
non-condensable gas at the maximal bubble radius, 7 is the adiabatic index of the non-condensable 
gas, p is the liquid density, and c is the speed of sound in the liquid. 



I. INTRODUCTION 

Research on cavitation is currently experiencing a re- 
birth within hydrodynamics. While traditionally associ- 
ated with turbine erosion [l| , cavitation bubbles are now 
exploited as tools in surgery [H, microchip cleaning 
water treatment [4] , and microfluidics [El, ^ . This wide 
spectrum of new applications relies on the diversity of 
processes associated with the collapse of cavitation bub- 
bles. Detailed studies revealed that these processes in- 
clude (i) the formation of rebound bubbles [3 , (ii) acous- 
tic shocks [H, (iii) thermal effects, leading to chemical re- 
actions [il and luminescence [lol^lili, and (iv) micro-jets 
[l2|-[l3- However, today there is no theory predicting 
the fraction of energy transferred into each of these pro- 
cesses. In the quest for such a theory, it seems wise to 
start with spherically collapsing bubbles, which produce 
no jets [14]. We also note that thermal processes typi- 
cally absorb negligible energy fractions 0. The problem 
then reduces to how the energy is split between rebound 
and shock in the spherical collapse. 

This paper presents an experimental and theoretical 
investigation of the energy partition between rebound 
bubbles and shocks. We first describe the experiment, 
which uses a mirror-focused laser pulse in combination 
with micro-gravity conditions to produce bubbles of ex- 
tremely high spherical symmetry. We then analyze mea- 
surements of the rebound sizes and the shock pressures of 
spherical bubbles produced in various experimental con- 
ditions. Interestingly, the energy ratio between rebound 
bubble and shock wave is found to decrease with the liq- 
uid pressure. We show that these observations can be 
explained using the Keller-Miksis collapse equation for 
a compressible liquid [l5| in combination with an adia- 
batic treatment of the non-condensable gas. Finally, we 
use this model to predict the energy partition between 



rebound and shock in a wide range of experimental con- 
ditions. 



II. EXPERIMENTAL SETUP 

The cavitation bubbles are obtained by focusing a 
high-energy laser in water (Fig. [1]). The laser source 
is a Q-switched Nd:YAG laser (Quantel CFR 400) deliv- 
ering pulses having 230 mJ maximum energy, 8 ns du- 
ration and a wavelength of 532 nm. The laser beam 
of 5 mm in diameter is expanded ten times before be- 
ing focused by an off-axis parabolic mirror with a fo- 
cal length of 54.5 mm, which is fixed inside the water 
container. The use of a parabolic mirror rather than 
an optical lens improves the focus by avoiding refraction 
and spherical aberration. The convergence angle is 53° 
and the focal point is located at the center of the water 
container. Owing to this large angle, the plasma gener- 
ated at water breakdown is more compact and spherical 
than in previous studies [l6|. The ensuing bubbles are 
so spherical that the faint hydrostatic pressure gradient 
due to gravity induces visible jets against the gravity vec- 
tor [1^. To avoid this source of asymmetry the experi- 
ment is performed in micro-gravity conditions [European 
Space Agency (ESA), 52nd parabolic flight campaign]. 
The flights consist of a total of 93 parabolic trajectories, 
flown by the aircraft A300 zero-g. Each parabola offers 
20 s of microgravity (acceleration < 0.01^). Given those 
unique conditions, our experiment produces millimetric 
cavitation bubbles of extremely high sphericity. 

The bubble is observed with a high-speed camera 
(Photron Fastcam SAl.l) at a rate of up to 250,000 
frames/s with an exposure time of 370 ns. A 3W light- 
emitting diode (LED) light with a small opening angle 
of 6° is used to illuminate the bubble from the back and 
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visualize the shock wave by shadowgraphy. The shock 
waves emitted at the generation and the cohapse of the 
bubble are monitored by a piezo-resistive dynamic pres- 
sure sensor. The pressure in the vessel is controlled by a 
vacuum pump. 




FIG. 1. (Color online) Sketch of the experimental setup 
aboard the A300 zero-g aircraft. 



III. EXPERIMENTAL RESULTS 

In the course of the flights, the bubble dynamics at 
three distinct water pressures Poo (10, 30, and 80 kPa) is 
observed. For each pressure, the laser pulse energy is var- 
ied from 55 to 230 mJ, resulting in maximal bubble radii 
Rmax from 2 to 5.6 mm. Figure [2] shows the normalized 
radius R/Rmax for a representative selection of bubbles 
as a function of the normalized time ^/tc, where Tc is 
the bubble collapse t ime ac cording to Rayleigh theory 
Tc = 0.915Rmax\/ p/ with p being the density of 
the liquid and Ap being the "driving pressure," i.e., the 
difference between the static liquid pressure p^o and the 
pressure py of the condensable vapor inside the bubble. 
The value of py is calculated with the Antoine equation 
from the temperature of the water measured for each of 
the three flight days. The three temperatures are, re- 
spectively, 16.8, 23.9, and 20.9 °C, corresponding to pv 
of 1910, 2950, and 2460 Pa. Ah the curves are remark- 
ably superposed during the first collapse, and are consis- 
tent with the Rayleigh theory (solid line in the figure). 
However, the dynamics of the rebound is very different 
depending on the pressure in the liquid p^o (see also Fig. 
[3|). The high-speed movies reveal that the normalized 
maximum radius of the first rebound bubble Rreb / Rmax 
decreases with poo- To interpret this result in terms of 
energy, we calculate the potential energy of a bubble as 

m 

Epot = / WApdr = —R^Ap, (1) 
Jo 



where R is the bubble radius. In particular, we define the 
initial bubble energy Eq and the rebound energy Ereb as 

47r 47r 

E0 = Y^max^P and Ereb = Y^reb^P' (2) 

Figure m reveals that Ereb scales with Eq for given liq- 
uid pressures Poo- 

To complete the picture, we consider the energy carried 
away by the spherical shock produced at the first bubble 
collapse. Given a shock pressure p(t), measured at a 
distance d from the bubble center, the shock energy is 
given by [l6| 

Esw = ^ / P{t)^dt, (3) 

where p is the water density and c is the speed of sound 
in water. In our experiment, the duration of the shock 
transition, i.e., the characteristic time scale ofp{t) (<100 
ns), is much shorter than the characteristic response time 
(10 ps) of the pressure sensor. Nevertheless, a rough es- 
timation of the shock energy remains possible under the 
assumption of a linear response. Explicitly, if we define 
h{t) as the sensor's impulse response, the response of the 
sensor s{t) is expressed as s{t) = h{t)^p{t)^ where de- 
notes the convolution. We assume that the pressure p{t) 
has a universal shape in the sense that pit) = PmaxP{t), 
where p{t) is the same function for all bubbles. The sig- 
nal can then be expressed as s{t) = h{t) ^PmaxP{t), and 
hence / s{t)dt = Pmax J h{t) * p{t)dt oc Pmax- In other 
words, Pmax IS proportioual to the integrated response. 
Substituting into Eq. (|3]), we finally obtain fil\ 

Esw oc J pl^^^p{tfdto<ipl^^^ (x(^j s{t)d-t^ . (4) 

The constant of proportionality in Eq. (|4]), which is 
unknown, is estimated such that the shock energy Esw 
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FIG. 2. (Golor online) The normalized radius for a repre- 
sentative selection of bubbles as a function of the normalized 
time, for different pressure levels poo- The experimental data 
(dots) are consistent with the Rayleigh theory (solid black 
line) . 
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FIG. 3. Selected high-speed images of a cavitation bubble at two different water pressures. The images are scaled so that 
the bubble appears with the same normalized Rmax on the figure, (a) poo = 30 Pa, Rreb/Rmax = 0.22 , (b) poo — 10 Pa, 

Rreb/R max — 0.57. 



equals the initial energy Eq in the extreme cases, where 
only a negligible rebound bubble is observed. The mea- 
sured shock energy Esw versus the initial potential en- 
ergy £^0 are presented in Fig. [5l Unlike the rebound en- 
ergy Ereb (Fig. H]), Esw is roughly proportional to Eq 
with no clear dependence on pressure. As we will show 
below, this is well explained by the fact that our exper- 
iments all lie in a "shock-saturated" regime, where the 
shock absorbs most of the available energy {Esw ^ Eq). 



IV. THEORETICAL MODEL 

Hereafter, a theoretical model is developed to com- 
pute the energies of the rebound bubble and the shock 
wave as a function of various experimental parameters. 
The standard model for the evolution of spherical cavi- 
tation bubbles, i.e., the Rayleigh-Plesset equation, can- 
not produce rebound bubbles and shock waves. To cal- 
culate the rebound motion it is necessary to include a 
non-condensable gas inside the bubble. We here assume 
that this gas is compressed and decompressed adiabati- 
cally, i.e., that is without heat transfer across the bubble 
surface. According to the adiabatic theory, the pressure 



Pg{t) of this non-condensable gas is then given by [18] 



Rr) 



37 



(5) 



where Pgo is the pressure at the maximal initial bubble 
radius Rmax-, R{t) is the evolving bubble radius, and 7 
is the adiabatic index also known as "heat capacity ra- 
tio." Second, to incorporate shock waves, we require a 
model for the bubble evolution in a compressible liquid. 
We here use the Keller-Miksis model [l5| , which is an ex- 
tension of the Rayleigh equation to compressible liquids, 
accurate to first order in the speed of sound c. As shown 
by Prosperetti [19] this model belongs to a more general 
class of first-order models and can be rewritten as 



R 



{pg - Ap){l + i;) + Rpg/c - (3 - v)R^p/2 
(1 - v)Rp 



(6) 



where v{t) = R{t)/c. Note that we deliberately neglect 
the effects of surface tension and viscosity for two rea- 
sons. First, these effects are quite irrelevant for the 
large bubbles in our experiment. Second, surface ten- 
sion and viscosity are generally insignificant at the last 
stage of the bubble collapse, since inertial forces increase 
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FIG. 4. Measured potential energy of the rebound bubble as 
a function of the initial bubble energy, for different pressure 
levels Poo . 



FIG. 5. Estimated energy in the shock wave as a function of 
the initial bubble energy, for different pressure levels poo • 
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more rapidly than viscous forces and surface tension as 
R{t) 0. The latter can therefore be neglected to cal- 
culate rebounds and shocks. 

Equations. (|5j) and (|6j), fitted with the initial conditions 
R{0) = Rmax. R{0) = 0, pg{0) = pgo, and pg{0) = 0, 
constitute a model for the collapse and the rebound of 
a spherical bubble, while including compression waves 
(shocks). We use the Runge-Kutta method to solve this 
model numerically. The radius R{t) is calculated as the 
bubble first collapses and then rebounds until it reaches 
its maximal rebound radius Rreb- 

Given a time solution of Eqs. (|5|) and (|6|) we can 
then calculate various energies. The initial bubble en- 
ergy £^0 and the energy of the rebound bubble Ereb are 
computed directly using Eq. (|2]). It is important to 
note that the temperature of the non-condensable gas 
changes during the adiabatic compression and decom- 
pression. The gas temperature at the rebound point is 
different from the initial temperature. Hence the internal 
energy U = {4:7r/3)R^Pg/{j — 1) of the non-condensable 
gas changes. We can calculate this energy change AU 
simply by subtracting the final value of U from the ini- 
tial one, 



AU = 



47r 



3(7-1) 



{PgoR 



3 

max 



Pg,rebRreb) 



(7) 



The adiabatic nature of the process implies that AU must 
be equal to the total work done by the liquid onto the 
non-condensable gas. This work can be calculated as 



AU = J SW = - J PgdV = - J AnR^Rpg dt, 



(8) 



where the time integral runs from the initial bubble ra- 
dius through the collapse point to the maximal rebound 
radius. To check the accuracy of our numerical solution 
we compute AU using both Eqs. ([7|) and (|8]). 

Given At/, the initial energy Eq^ and the potential en- 
ergy of the rebound E^eb^ the compression energy of the 
shock wave Esw can be computed from energy conser- 
vation as 



AU 



(9) 



Esw = Eq — Ereb 

Finally, we introduce the energy fractions 

Creb = Ereb/Eo, esw = Esw/Eo, eu = AU / Eq. (10) 

Equation ([9]) implies the normalization ereb^^sw = 
1. 

How do Creb-, ^sw-, and eu depend on the six model 
parameters Rmax-, Ap, pc/O, 7, and c ? We first note 
that the four energies Eq^ Ereb^ Esw^ and AU all scale 
as Rfnax- This can be shown by rewriting the model as 
a function of the normalized radius r{t) = R{t) / Rmax- 
Therefore Crebi ^swi and eu are independent of Rmax- To 
test the remaining five model parameters we ran 2.7 x 10^ 
independent computations of Crebi ^sw-, and eu by tak- 
ing logarithmically spaced parameters from the following 



intervals: Ap e [1,100] kPa, Pgo e [0.1,100] Pa, p e 
[500,1500] kgm-^ c e [1000,2000] ms-\ 7 G [1.3,1.5]. 
By systematically studying the variation of ereb, ^sw, 
and eu as a function of the five parameters, we can draw 
two main conclusions. First, the internal energy fraction 
is negligible because eu < 0.01 in all situations. Second, 
all variations of ereb and esw as a function of the five 
model parameters Ap, pgo, 7, p, and c can be explained 
by using a single non-dimensional parameter 



1-1/7 ■ 



(11) 



In fact. Fig. (j6j) shows the 2.7-10^ values of ereb and esw 
as a function of ^, revealing a tight correlation. The 
parameter ^ was found by first constructing the non- 
dimensional parameter Ap {pgo)^{pc^)~^~^ from the four 
dimensional parameters Ap, pgo^ p, and c. The computed 
results are then grouped depending on the value of 7. a 
is determined for each group as the value that maximizes 
the Pearson correlation coefficient for ereb ^ [0.2,0.8]. 
We restrict ereb to the interval where a small variation in 
^ leads to a large variation of e^eb, thus where we want 
the relation to be the most univocal. The values of a 
obtained depend on 7 as a = I/7 with an error of ±10%. 
As the curves obtained for each value of 7 are still hori- 
zontally shifted, the second step is to introduce the factor 
7^ . /3 = 6 is then determined by maximizing the Pearson 
correlation coefficient on ereb ^ [0.2,0.8] for all values of 
7- 
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FIG. 6. (Color online) Fraction of energy in the rebound 
Creb and in the shock wave esw as a function of the non- 
dimensional parameter ^. The solid curves are the results 
from the theoretical model. The discrete black symbols are 
the values obtained experimentally, along with the measure- 
ment error bars. The white symbols are data extracted from 
the literature. 
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V. DISCUSSION 

A. Comparison between model predictions and 
experiment 

The theoretical model allows us to explain why, ac- 
cording to our experimental results, the energy of the 
rebound depends on the pressure in the hquid while the 
energy of the shock wave seems to scale with the initial 
potential energy only. However, to compare the experi- 
mental results with the theoretical ones, we need a value 
for pgo in addition to the measured Ap and Rmax and 
the known p, c, and 7. Since Pgo is not directly mea- 
surable, we simply assume this pressure to be constant. 
We estimated its value by fitting the model to the ex- 
perimental results. For each measurement, the value of 
PgO leading to the observed Rreb is calculated with an it- 
erative process. The results are averaged and we obtain 
p^ = 7.0 ± 3.5Pa. The relatively small variance a poste- 
riori justifies the assumption of a constant value for pgo. 
The experimental points are plotted in Fig. [6l where the 
values of ^ are calculated using p^. We observe that all 
our experimental data lies in a regime where esw ^ 1- 
So when ^ varies because of the change of Ap, the rela- 
tive difference is important for Creb but not for esw- As 
Creb and esw represent the slopes of the curves in Figs, 
m and m respectively, the difference in poo is significant 
for the rebound, but insignificant for the shock. 

We observe, both theoretically and experimentally, the 
relation Ereb^Esw = ^o- The results using the theoreti- 
cal model show that /SXJ in Eq.(j9]) is negligible, which im- 
plies Ereb^Esw = Eq. And when adding the experimen- 
tal data from Figs. |4]and[5l we obtain Ereb-^ Esw = £^0, 
within a scatter of ±20%. This scatter is reasonable con- 
sidering the uncertainty introduced with the rough esti- 
mation of Esw- 



B. Comparison with earlier work 

The main issue with the treatment of the rebound is 
that the pressure of non-condensable gas pgo is usually 
needed in the equation of motion. As pgo is not mea- 
surable and its origin is not clear yet, it is difficult to 
estimate it and thus to validate a model. So in a concern 
of evaluating our theoretical model, we look at previous 
studies for comparison. In the experimental and numer- 
ical work of Kroninger et al. |2Q|, and in the numerical 
work of Sadighi-Bonabi et al [21| and Fujikawa et al 
^2], we found estimates of pgo or enough information to 
obtain them. The data extracted from these articles are 
plotted in Fig. [6l We observe that despite the differ- 
ent treatments of the thermodynamics inside the bubble 
(Ref. [20I used a van der Waals equation, Ref. [Ill con- 
sidered hydrochemical reactions, and Ref. [22| considered 
conductive heat transfer and condensation or evapora- 
tion) our model reproduces reasonably well their results. 
Yet, the drawback of our model is that the temperature 



at the collapse is overestimated because of the neglected 
thermal transport. This could be improved by the ad- 
dition of heat transfer or chemical reactions, but at the 
cost of the simplicity of the model. 

Akhatov et al propose a mathematical model sup- 
ported by experimental measurements of the rebound of 
a spherical cavitation bubble. Because of the difference 
in the model used (the pressure of non-condensable gas 
is derived from a phase transition equation), we could 
not derive a value for pgo for quantitative comparison of 
our results. Nonetheless, qualitatively, the conclusions 
are the same. Akhatov et al. observed that the ratio 
between the radius of the rebound and the initial bubble 
is constant when only varying the initial radius of the 
bubble, which confirms the univocal relation between ^ 
and Creb' They also showed numerically that for given 
experimental conditions, when the concentration of the 
non-condensable gas in the bubble increases, the radius 
of the rebound bubble increases too. This also agrees 
with our conclusions. Indeed, the increase of the concen- 
tration of non-condensable gas means a smaller value of 
^ which implies, according to Fig. [6l an increase of €reb 
and thus of the rebound radius. 



C. Negligible role of gravity 

We have already demonstrated (see Obreschkow et al. 
Q) that gravity can affect the collapse of a cavitation 
bubble in the form of the occurrence of a vapor jet (see 
Fig. [7]). The volume of the vapor jet normalized to the 
maximum volume of the rebound was found to be propor- 
tional to the non-dimensional parameter ( = \\/p\Ro/Ap, 
where Vp is the hydrostatic pressure gradient. Therefore 
we also performed the experiments in this paper with the 
same parameters at normal gravity (1^) and hypergrav- 
ity (1.8^). The values of the non-dimensional parameter 
were C e [2.5 • IQ-^, 7 • 10-^]. Unlike for the vapor jet, we 
do not observe a significant difference in the energy par- 
tition when the gravity changes. The relative difference 
between the values of Creb at 0^ and at >1^ are smaller 
than the standard deviations of the measurements at 0^. 
We deduce that the energy transferred into the vapor jet 
is negligible compared to the energy in the rebound and 
in the shock. In consequence, the results of this paper 
also apply to bubbles collapsing in a hydrostatic pressure 
gradient for < 7 • 10~^. 

Note that in the studies cited here 0, Hoj-lil and 
plotted in FiglHl the non-dimensional parameter is C ^ 
7 • 10-^: We found respectively 1.5 • IQ-"^, 7.5 • 10-^ 
3.0 • 10~^, and 1.4 • 10~^. We thus consider that the com- 
parison of our results with these data is justified. 



D. Implications 

A systematic experimental and theoretical investiga- 
tion of the rebound and shock energy at the collapse of 
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FIG. 7. Collapse and rebound of a bubble in 0^ (upper) and 
Ig (lower). Note the shock visible at the collapse, and the 
vapor jet on the rebound for Ig [3]. 



a spherical cavitation bubble is presented. This led us 
to identify a single non-dimensional parameter ^, which 
links the experimental conditions to the fraction of en- 
ergy in the rebound bubble and in the shock wave. This 
finding has important implications for many engineering 
applications. Depending on the desired requirements, we 
can tune the experimental parameters in order to obtain 
a value for ^ that favors either the rebound or the shock. 
For example in micro-pumping, an enhanced rebound is 
desired to increase efficiency, because the volume of the 
bubble affects the displacement of the liquid [6[. Con- 
versely, for applications that rely on cavitation erosion, 
stronger shocks would be preferred to accelerate the pro- 
cess. 

The methodology presented here to estimate pgo also 
opens perspectives for understanding the origin of the 
non-condensable gas in the bubble at the collapse. So far 
the non-condensable gas has been assumed to be a combi- 
nation of trapped vapor, laser breakdown products, and 
gas initially present in the water [7]. A method to ver- 
ify this would be to systematically vary the experimental 
conditions and assess their effect on pgQ. Our model could 
then be used to extract the values of Pgo by fitting the ex- 
perimental results with the theoretical ones. In the same 
line of thought, we could estimate water properties, such 
as the concentration of dissolved gas and nuclei, solely 
based on observing the rebound of natural or artificially 
generated cavitation bubbles. These results could be 
complemented with observations of extreme cases such 
as cryogenic fluids, where the experimental parameters 
("driving pressures" Ap, adiabatic index 7, density p. 



and speed of sound c) are very different compared to the 
case of a bubble in water [23, S3 • 



VI. CONCLUSION 

We observe experimentally that the pressure of the liq- 
uid affects the energy partition between the rebound and 
the shock wave at the collapse of a highly spherical laser- 
induced bubble. This unique experiment is performed in 
microgravity conditions in order to avoid the formation 
of a microjet due to the hydrostatic pressure gradient 
at the collapse of the millimetric bubble. A theoretical 
model for the collapse of spherical bubbles is proposed, 
relying on a compressible equation of motion and an adi- 
abatic treatment of the non-condensable gas inside the 
bubble. The partition of the energy between rebound 
and shock is calculated for a wide range of parameters. 
It appears that, in addition to the pressure in the liq- 
uid, the physical properties of the liquid and the pressure 
of non-condensable gas inside the bubble affect the en- 
ergy partition. These parameters can be combined into 
a single non-dimensional parameter ^, which dictates the 
energy partition. 

The ability to predict the energy partition between re- 
bound and shock is valuable in many engineering applica- 
tions. The operating conditions can be adjusted to favor 
rebound or shock depending on the requirements. Using 
the method developed in this paper, it is also possible to 
estimate the pressure of the non-condensable gas in the 
bubble by fitting, with an iterative process, the experi- 
mental observations of the radius of the rebound bubble 
with the theoretical results. However, as the pressure 
of the non-condensable gas is not directly measurable, 
this method still has to be quantitatively validated with 
further experiments. 
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